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© On-line determination of polymer properties in a continuous polymerization reactor. 



© The molecular weights (Mw and Mn) and polydispersity (Mw/Mn) of a polymer have significant impact on the 
final product properties and processability of materials made out of it. This invention provides a method for on- 
line prediction of polymer properties. In a continuous polymerization reactor in which an irreversible polymeriza- 
tion reaction occurs, sensors provide axial temperature and pressure profile in the reactor as well as monomer 
flow rates and concentration for determining Mw, Mn and polydispersity. Based on a non-Newtonian fluid 
mechanics analysis, differential influx analysis of a side stream incorporation and determination of the moments 
of the molecular weight distribution via a Newton-Raphson algorithm, the model calculates an axial polymer 
properties profile. 

On-line real time "analytical measurement" of polymer properties is used to alert the operator of unusual 
reactor problems, to optimize reactor operation during shut-down and start-up or rate changes. This allows 
correlation of polymer properties to product performance in downstream processes. 
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On-Line Determination of Polymer Properties In a Continuous Polymerization Reactor 



Background of the Invention 

A major problem constantly faced by engineers in the polymer reaction engineering field is to 

5 accurately predict the performance of continuous polymerization reactors. Due to the critical commercial 
application of polymers, the polymerization reaction should be very carefully controlled in order to insure 
uniform properties. However, most commercial installations lack proper on-line techniques to predict the 
product behavior as it is made in the reactor. A number of parameters such as hold-up time, pressure and 
temperature profile in the reactor are key factors affecting reactor and product performance. Also, product 

70 properties are affected by reactor start-up, shut-down and rate changes. There is indeed the need for a 
quick, continuous on-line measurement of polymer properties. As an on-line analyzer of polymer properties 
has not been demonstrated, the approach of using computational fluid dynamics techniques, elementary 
measurement of temperature, pressure and flow by means of in-line instruments and computer hardware 
(DCS, VAX 8000) have proven successful. Off-line verification of model predictions by means of a Liquid 

75 Chromatograph is then easily done. 

This invention provides on-line real time polymer product properties as the polymer is being made in 
the reactor without making use of the chemistry of the reacting system. In fact to follow the kinetics of the 
reaction for making polymer properties prediction has proven to be a approach without recognizeable value, 
mostly for the reason that energies of activation and reaction rate constants for the reaction system are hard 

20 to measure accurately and, therefore, the numerical solution of the set of differential equations is inaccurate. 
Also, the polymerization chemistry might not be well understood, in particular, the kinetics of cyclic dimer 
and oligomer formation. 

The present invention deals with the application of fundamental principles of fluid mechanics and 
polymer processing that were used in the on-line prediction of molecular weight weight average (Mw), 
25 molecular weight number average (Mn) and polydispersity (Mw/Mn) in a continuous tubular polymerization 
reactor. 



Summary of the Invention 

30 

The main objective of the present invention is to provide a method for monitoring and controlling the 
molecular weight build-up of a polymer inside a continuous polymerization reactor. The number average 
and weight average molecular weights are predicted on-line from fluid dynamic principles and a slow influx 
analysis of multiple side streams incorporation into the main bulk of the fluid flow. Process sensors placed 

35 along the length of the reactor provide continuous on-line measurement of the input variables required to 
perform the model calculations. 

In a given reaction section, information regarding the average temperature, pressure drop per unit of 
length, fluid shear rate and residence time are coupled with the mixing times required to obtain a given side 
stream incorporation to determine polymerization solution viscosities. The solvent viscosity and solution 

40 viscosity determine a specific viscosity that in conjunction with information on polymer fraction in solution, 
solution density and polymer concentration are used to determine the polymer solution "intrinsic viscosity". 
In this step in the solution, a Newton-Raphson convergence algorithm is required. The molecular weight Mw 
of the polymer in the particular reactor section is obtained from the value of intrinsic viscosity via the Mark- 
Houwik equation. 

45 The determination of the molecular weight number average Mn can proceed by first calculating an 
external stoichiometric ratio from values of the monomer's flow rates. Subsequently, the first and second 
moments of the molecular weights of the individual monomers are obtained. Then, a guess on the apparent 
stoichiometric ratio is made on the basis of the external stoichiometric ratio. Next, the calculation of the first 
and second moments of the molecular weight distribution of the polymer can proceed from knowledge of 

so the apparent stoichiometric ratio and the first and second moments of the. individual monomers. At this 
point, the molecular weight weight average Mw is calculated as the ratio of the second moment to the first 
moment of the molecular weight distribution of the polymer. If the value of Mw calculated in this fashion is 
the same as the value of Mw calculated in the paragraph above, then the solution has converged and we 
have at hand the right numerical values of the moments of the molecular weight distribution of the polymer. 
Therefore, the value of molecular weight number average Mn is obtained from the first moment of the 
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molecular weight distrioution of the polymer and the apparent stoichiometric ratio. In the event that 
convergence was not achieved, a Newton-Raphson solution guesses a "new value" of the apparent 
stoichiometric ratio and repeats the logic explained in this paragraph until convergence Is achieved. 

The determination of the polydispersity of the* polymer can proceed in a straightforward fashion by 
5 using the converged values of the first and second moments of the molecular weight distribution of the 
polymer. The polydispersity of the polymer in a given section of the reactor is equal to the ratio of the 
second moment of the molecular weight distribution to the square of the first moment of the molecular 
weight distribution. 

The present invention provides real time information on the reactor performance and makes it possible 
jo to implement a closed loop control of molecular weight number average Mn and polydispersity. Since the 
reactor control via an on-line viscometer provides means only to somewhat control molecular weight weight 
average Mw (see Fig. 4), the approach developed in this invention would without a doubt provide means for 
. improved control of a continuous polymerization reactor. In fact, viscosity might not be a suitable control 
variable. In certain systems, viscosity may remain essentially unchanged while the molecular weight 
is number average Mn varies out of control (see Fig. 5). This invention becomes an on-line real time 
"analytical measurement" tool of polymer characteristics. In fact, off-line laboratory analysis of polymer 
solution samples taken along sample ports in the reactor will take too long (more than 24 hours) versus on- 
line prediction. The model can make predictions as frequently as desired (as low as every 5 seconds) in as 
many reactor sections as appropriate input data is available. This invention allows the on-line monitoring of 
20 product quality in the continuous reactor while avoiding time-consuming laboratory analysis. The calculated 
results can then be used to alert the operator of unusual reactor problems, to optimize reactor operation 
during shut-down, and start-up or rate changes: Since all of the calculated results are stored in the data 
acquisition system, this invention allows the future correlation of product performance in downstream 
processes to the polymer characteristics. 

25 

Brief Description of the Drawings 

Fig. 1 shows a flowsheet of the continuous mutti-pais polymerization reactor system. 
30 Fig. 1 A is similar to Fig. 1 showing a flow sheet of a three pais reactor system. 

Fig. 2 shows a block diagram of the computer hardware and logical process information flow. 
Fig. 3 shows the flowchart explaining the logical software calculation steps. 

Fig. 4 shows typical experimental values of on-line control viscometer versus molecular weight weight 
average Mw. 

35 Fig. 5 shows typical experimental values of on-line control viscometer versus molecular weight 

number average Mn. 

Fig. 6, 7 and 8 are graphs of comparing model prediction with experimental results obtained in 
Example II. 

Fig. 9 is a graph comparing bulk viscosity with model prediction in Example I. 

40 

Detailed Description of the Preferred Embodiment 

The following contains a thorough description of the invention with regards to the on-line prediction and 
control of the polymer characteristics in a continuous polymerization reactor in which an irreversible 

45 polymerization reaction occurs. 

Fig. 1 depicts the flowsheet regarding the present invention. A monomer "B" dissolved in' a solvent 
contained in storage tank T1 is pumped by means of a pump P1 through a flow-metering system M1 into 
the inlet port of a continuous polymerization reactor R1 . At the same time, pure monomer "A" contained in 
storage tank T2 is pumped by means of pump P2 through a series of fiowmetering systems "mi" that are 

so located at different axial locations in the reactor. The individual fiowmetering systems "mi" allows the 
addition of varying amounts of monomer "A" into the reactor to optimize temperature profile and facilitate 
thermal management of the exothermic or endothermic reactions. As monomer "B" flows inside the reactor, 
it reacts with side-streams that inject monomer "A". The over-all scheme of the polymerization can be 
represented by: 

55 [B-B] + r[A-A] > [AA-BB]n-AA + Byproducts (1) 

where r < 1.0 and n = 1, 2, 3 

r [ = ] overall stoichiometric ratio. 

At the exit of the continuous polymerization reactor, the polymer and byproducts are stored in storage 
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tank T3 and pump when necessary for further chemical processing using pump P3. A coolant stream S9 is 
fed to the reactor and exits at S10 for coolant or heating purposes in either a cocurrent or countercurrent 
fashion. The flowrates of monomers and coolant are controlled appropriately by means of flowmeters to 
provide desired stoichiometric ratios or to control the reactor within a given temperature profile. 

5 Rg. 1A illustrates schematically a three pais reactor R2 utilizing the same general flow system as 
described above. The coolant flow has been deleted. The invention may be operable without a cooling 
jacket depending on the polymer. The pressure and temperature sensors designated T1-T6 and P1-P6, 
respectively, are shown at the inlet and outlet ends of each reactor pass, i.e. pass 1 , pass 2 and pass 3. 
The analyzer CB-1 for monomer B concentration is located near the outlet end of pump P1. The tachometer 

w sensors RPM-1, RPM-2 and RPM-3 are associated with motors EM-1, EM-2 and EM-3, respectively. The 
motors drive the agitators positioned in the respective passes of the reactor. 

A list of instrumentation used is presented in the following Table 1. AJJ of the sensors, flow metering 
devices, control valves and concentration analyzers are commercially available items and are used at 
locations indicated in Figs. 1 and 1A. 

TABLE 1 





TYPICAL LIST OF INSTRUMENTATION USED 


20 


ITEM 


MODEL 


COMMERCIAL 
IDENTIFICATION 




SENSORS 


25 
30 


RTD Sensors (T1, 
T2, T3. T4, T5. T6) 
Pressure Sensors 
(P1.P2, P3, P4, P5, 
P6) 

Tach Sensors 

(RPM1.RPM2, 

RPM3) 

Cone. Analyzer 
Monomer B In 
Solvent CB 


RTD-18-HT-9-OT-1 Sheath 

Material Tantalum 

YSTG1 40E1 600000-MBFIC3 

Part No. ARCJ-2A-001 A and 
MP-37CA 

Model 1801 


Thermo Electric Saddlebrook, 
NJ 

Honeywell Fort Washington. 
PA 

Red Lion Controls Willow [ 
Springs Circle York, PA 

E. 1. du Pont de Nemours and 
Company Wilmington. DE 


35 


FLOW CONTROL 




Turbine Meters 


40 
45 


Monomer B M1 
Monomer A M2 
Monomer A M3 
Monomer A M4 


FT-6-16-CX50-UC 
FT-6-8CX5-UC 
FT0-CX10-UC 
FT0-CX6-LJC 


Flow Technology, Inc. 

Phoenix, AZ 

Flow Technology, Inc. 

Phoenix, AZ 

Row Technology, Inc. 

Phoenix, AZ 

Row Technology, Inc. 

Phoenix, AZ 




Control Valves 


50 
55 


Monomer B Ml 

Monomer A M2 

Monomer A M3 
Monomer A M4 


1456-2 Size 1" Trim. 3/18" Un 

785-73B Size 1/2" Trim P 

785-73B Size 1/2" Trim I 
785-73B Size 1/2" Trim K 


Kieley & Mueller Middletown, 
NY 

Badger Meter, Inc. Research 
Control Valves Tulsa, OK 

Badger Meter. Inc. Research 
Control Valves Tulsa. OK 



Rg. 2 shows the block diagram of the computer hardware and logical process information flow. Each 
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process sensor 10 (Temperature TS, Pressure PS, Row FS, Tachometer RPMS and Monomer Concen- 
ration CM) is connected to an analog to digital converter card 20 in the input/output (I/O) nest of the unit 
control module 30 (UCM) of the distributed control system 40 (DCS). When this hook-up is done and the 
cards are properly calibrated, the UCM must be configured to read the digital values 20 into its memory 60 

s and convert them to the desired range for calculation and display 50 at the DCS video spec operator's 
console. The live process input data required by the model must next be transferred to the VAX 8000 
computer 100. This is done by the Foxboro computer interface 80 (FCI) which communicates with the 
technical development system 90 TDS. A VAX procedure runs the model subroutines 130 at a specified 
frequency set-up by the time scheduler 120. The time scheduler also hibernates the model subroutines. 

w The input data required by the model subroutines is provided by the Data Acquisition System 110. The 
model subroutines return the calculated values of polymer properties in a reactor section back into the Data 
Acquisition System where loops are assigned in the Data Acquisition System for display 140 or to be kept 
in memory 150. The calculated results can be loaded down from the Data Acquisition System 110 through 
the FOXBORO computer interface 80 back into the memory 60 or directly into a controller card 70 in the 

75 UCM for further download into the control valves 160. 

Fig. 3 presents the flowchart explaining the logical software calculation steps in the execution of the 
model subroutine. After the time schedule 120 has called the wake-up subroutine from hibernation FC1, the 
model subroutine obtains the input parameters FC2 (pressure, temperature and flow of monomers) in a 
given reactor section and checks for non-zero values of input parameters FC3. If there are indeed zero 

20 values, an alarm FC4 is displayed and the subroutine goes back to FC2 and stops calculations. If the input 
parameters are fine, then an average temperature, pressure drop and residence time in the reactor section 
are calculated FC5. The procedure continues on to determine the mixing times required in the fractional 
incorporation of a side stream of monomer "A" into the main bulk of fluid FC6. 
F = 1.-EXP(al*(Trrb!) (2) 

25 where 

al = F(% incorporation of monomer "A") 
bl = F(al, flow rate) 

Tr = residence time from inlet to section in reactor 

The solvent viscosity in the reactor FC7 is determined by means of the following equation: 
30 Usot = cl*exp(dl/(el + t) (3) 
where 

cl = constant dependent on solvent type 

dl « temperature dependent constant 

el = temperature dependent constant 
35 t = average temperature in reactor section 

The solution viscosity in the reactor can be calculated from fundamental principles of Non-Newtonian 

fluid mechanics FC8. The key determination lies in a good prediction of the shear stress for flow inside a 

pipe reactor which can contain or not a shearing causing device like an agitator. This invention will cover 

either since the modification to the shear stress equation can be readily accomplished. The fundamental 
40 relationship for shear stress in a pipe is: 

Trz = (AP*D/4ir(r/R-(R/r)*X w 2) (4) 

where 

AP = pressure drop 

L = length of reactor section 
45 r = variable radius inside the flow field 

R = reactor inside radius 

X = constant of integration 

Since for Non-Newtonian fluids there are a number of models that relate the shear rate to the shear 

viscosity of the polymer, one can choose the one which is the most appropriate for the situation at hand. In 
so this invention we chose the Power Law model for Non-Newtonian fluids, represented by 

Trz = (u/gc)*(dv/drrn (5) 

where 

u = shear viscosity 
gc = gravitational conversion factor 
55 dv/dr = shear rate 

n = pseudoplasticity index 

The solution viscosity is obtained by setting equation (5) equal to equation (4) and then solving for the 
unknown u. In some instances, the shear viscosity in equation (5) can be made a function of the type of 
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agitator and the rotational speed of it. Also, the shear rate could be a function of the flow, rate only, such as 
flow inside a pipe without an agitator. 

The specific viscosity is then calculated by using the well-known definition in polymer chemistry 
according to FC9: 
5 Us - (u/Usol) (6) 

The calculation of polymer fraction in solution from monomers' flow rates is performed in FC10. Then, 
the calculation of solution density FC1 1 from average temperature in a reactor section and polymer fraction 
in solution is done. The equations herein mentioned will vary from one polymer system .to another and are 
proprietary. 

w The following step is to calculate the polymer concentration in solution FC12 from solution density and 
polymer fraction in solution. The equation is: 
C = 100*WTIO (7) 
where 

C = polymer concentration in solution, g/dl 
is W = polymer fraction in solution, % by weight 
RO = solution density, g/cc 

Calculation of the axial external stoichiometric ratio in every reactor section FC13 is accomplished from 
momomers' flows and the slow influx incorporation fraction. The external stoichiometric ratio is given as 
Re = sum(Fi monomer "A")/FB (8) 
. 20 where 

R = flow rate of monomer W A" injected in the i-th position along the reactor length 
FB = flow rate of monomer B B W fed to reactor inlet 

The prediction of the intrinsic viscosity [n] FC15 from the specific viscosity and pojymer concentration 
can proceed by a Newton-Raphson FC18 solution of the Martin equation: 
25 Us = C[nrEXP(K"-[nrC) (9) 

After convergence is achieved FC17, the cal culation of molecular weight weight average F19 can be 
obtained from the Mark-Houwik equation: 
Mw = (Us/ani-'b) (10) 
where 

30 a = constant dependent on polymer type 

b = constant dependent on polymer type 

The calculation of molecular weight number average Mn FC20 is begun by first assuming an apparent 

stoichiometric ratio from the external stoichiometric ratio based on ingredient flows in every reactor section. 

Ra=k'Re (11) 
35 where k< 1.0 

The next step is to determine the first and second moments of the molecular weight distribution of the 
individual monomers FC21 . If the monomers are pure (monodisperse). then the first and second moments 
of the molecular weight distribution of the monomers become equal to the actual molecular weight (based 
on summation of individual atomic weights of the atoms composing the molecule) of the monomer. 

40 The calculation of the first and second moments of the . molecular weight distribution of the reacting 
mixture FC22 in a reactor section is based on the apparent stoichiometric ratio and the first and second 
moments of the individual monomers. This calculation is performed by using the definition of the moments 
for a generalized polycondensation reaction: 
first moment 

45 M1 = (maa(1) + rmbb(1))/(1-r) (12) 
where 

maa(1) = first moment of individual monomer "A" 
mbb(1) = first moment of individual monomer "B n 
r = apparent stoichiometric ratio 
so second moment 

M2 = maa(2) + (r/(1 -r))-[{maa(1 )}"2 + 4*maa(1 )*mbb(1 ) + {maa(1 )}"2 + {mbb(1 )}~2] + [2Tqe/({1 -qe}*{1 -r})]*- 

[maa(1) + mbb(1)r2 (13) 

where 

maa(2) - second moment of individual monomer "A" 
55 mbb(2) = second moment of individual monomer n B w 
qe = effective stoichiometric ratio = r/Re 

At this point, the molecular weight weight average Mw is calculated as the ratio of the second moment 
to the first moment of the molecular weight distribution FC23 of the polymer. The equation is given as: 
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Mw = M2/M1 (14) 

If the value of Mw calculated by equation (14) Is the same as that calculated In equation (10) above, 
then the solution converged FC24 and we have at hand the right values of the moments of the molecular* 
weight distribution of the polymer. Therefore, the value of molecular weight number average FC27 of the 

s polymer is given by the first moment of the molecular weight distribution as given in equation (12). If, on the. 
other hand, convergence was. not achieved, a Newton- Raphson solution algorithm guesses a new value of 
the apparent stoichiometric ratio FC25 and repeats the calculation sequence explained in this paragraph 
until convergence is achieved. 

Determination of the polydispersity is. therefore, a straightforward calculation by using the converged 

w values of the first and second moment of the molecular weight distribution FC28. The polydispersity of the 
molecular weight distribution of the polymer is then: 
PF = M2/(M1)**2 (15) 

After final calculation of polydispersity. the model subroutine returns control to the scheduler FC29 
subroutine and can either go back and do the same calculations for another reactor section or stop 
is altogether and activate the hibernation subroutine FC30. Finally, the whole process will be repeated after 
the hibernation time limit is reached and then the model subroutine is called again and the calculation 
procedure is repeated all over. 

20 EXAMPLE I 

The continuous polymerization reactor system similar to that shown in Fig. 1A f except that five rather 
than three passes are used, is controlled by means of a closed loop system that senses the polymer 

25 viscosity at the exit of the fifth reactor pass and corrects the specified viscosity value by adjusting the flow 
of monomer n A w to the middle of the third pass. 

According to the input data required by the on-line math model, inlet and outlet pressure, temperature, 
monomer flows and rpm of the agitator device in the fifth pass were used to calculate the actual molecular 
weight weight average of the polymer for a three day period. The calculated Mw results appear in Fig. 9 in 

30 addition to the bulk viscosity (curve A) measured by the on-line viscometer at. the end of the fifth pass. The 
Mw is the average for the fifth pass as predicted by the model. It is worthwhile to notice that the general 
trend of the model Mw predictions follow the bulk viscosity because of their close relationship. A good 
example of this is the rightmost part of Fig. 9, just after a reactor start-up. However, in the short term, the 
model is predicting changes in Mw not shown by the bulk viscometer, hence suggesting perhaps that 

36 viscosity only might not be a suitable control variable. In certain systema, viscosity may remain essentially 
unchanged while the molecular weight number average Mn varies out of control. 

EXAMPLE II 

40 

\ 

This example deals with the evaluation of the polymer made during a two hour period in the continuous 
polymerization reaction shown in Fig. 9. Off-line verification of model predictions by analyzing polymer 
samples taken from the exit of. the fifth reactor pass was accomplished with a liquid chromatograph. 

4S Polymer samples were collected every three minutes, analyzed and compared to the on-line real time 
model predictions. Figs. 6, 7 and 8 show the molecular weight weight average, molecular weight number 
average and polydispersity, respectively. Table 2 shows the individual values of the polymer properties as 
well as the percent error of the model predictions versus the measured values. The data suggests good 
agreement between model predictions and measured values. In some instances, the differences could be 

so explained due to the fact that the sample port was located at the exit of the fifth reactor pass and the actual 
model predictions represent the average for the fifth pass. As proven herein, this method allows the on-line 
monitoring of product quality in the continuous reactor while avoiding time-consuming laboratory analysis. 
The calculated results can then be used to alert the operator of unusual reactor problems, to optimize 
reactor operation during shutdown and startup or rate changes as well as determine when to segregate 

55 product which is off-spec. Since all of the results are stored in the data acquisition system, this patent 
allows for the future correlation of product performance in downstrean processes to the polymer characteris- 
tics. 
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Claims 

1. A method of predicting physical properties of a polymer formed from a plurality of monomers flowing 
in a path through a continuous polymerization reactor in an irreversible reaction process with the aid of a 
digital computer, said method comprising: 

(a) obtaining a data set of flow rates of the monomers, pressure and temperature of the polymer at a 
plurality of locations along said path all on a real time basis; 
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(b) providing the computer with a data base that includes the data set obtained by step (a); 

(c) programming the computer according to a* Newton-Raphson algorithm; 

(d) calculating the physical property of the polymer based on non-Newtonian fluid mechanics 
according to the algorithm of step (c) using the data base of step (b); 

(e) displaying said property as a display element whereby a predicted physical property of the 
polymer is obtained based on non-Newtonian fluid mechanics. 

2. The method of claim 1 wherein said physical property is molecular weight. 
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IGURE 3. FLOW CHART EXPLAINING THE LOGICAL 
SOFTWARE CALCULATION STEPS 
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FIGURE « 4 

i 1 VISCOMETER READING VERSUS EXPERIMENTAL MOLECULAR WEIGHT Mw 
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FIGURE t S 

| 1 VISCOMETER READING VERSUS EXPERIMENTAL MOLECULAR WEIGHT Mn 
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